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The internal energy of hard spheres (HS) is the same as that of an ideal gas, so that the energy route to thermodynamics becomes 
useless. This problem can be avoided by taking an interaction potential that reduces to the HS one in certain limits. In this paper the 
i-^H ' square-shoulder (SS) potential characterized by a hard-core diameter a' , a soft-core diameter a > a' and a shoulder height e is considered. 

The SS potential becomes the HS one if (i) e — > 0, or (ii) e — > oo. or (iii) a' — * a or (iv) a' — > and e — » oo. The energy-route equation 
of state for the HS fluid is obtained in terms of the radial distribution function for the SS fluid by taking the limits (i) and (ii). This 
C , equation of state is shown to exhibit, in general, an artificial dependence on the diameter ratio a' /a. If furthermore the limit a 1 /a — » 1 is 
I , taken, the resulting equation of state for HS coincides with that obtained through the virial route. The necessary and sufficient condition 
, to get thermodynamic consistency between both routes for arbitrary a' /a is derived. 

CO 
, 

. 1 Introduction 

S: . 

i . As is well known, there exist several routes to obtain the thermodynamic quantities of a fluid in equilibrium 

£H I in terms of the pair interaction potential cj)(r) and the radial distribution function g(r) [1-3]. The most 

Q ■ frequently used are the virial route, 
O ' 



> 

■ the compressibility route 



^ - Z = 1 + \p /°°drr 3 y (r)|-e-^), (1) 
P "J jo 
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and the energy route 



1 2(3 
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n=A + 2vrp/ drr 2 <t){r)e-^y{r). (3) 



o 



In Eqs. P is the pressure, p is the number density, (3 = l/ksT is the inverse temperature, Z is 

the compressibility factor, x is the (dimensionless) isothermal compressibility, u is the internal energy per 
particle and y{r) = exp\ft(f)(r)]g(r) is the cavity (or background) function. 
The three thermodynamic quantities Z, \ and u are connected by thermodynamic relations, namely 

X" 1 = ^ ( P Z) , (4) 
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Thus, the compressibility factor, and hence the equation of state (EoS), can be obtained from y(r) [or, 
equivalently, g(r)] either directly from Eq. ©, or from Eqs. (J2J) and (J3J), or from Eqs. © and ©. Given an 
interaction potential (j)(r), if the exact cavity function y(r) is known for any thermodynamic state (p,/3), 
the three routes yield of course the same EoS. On the other hand, if an approximate function y(r) is used, 
a different result is in general obtained from each route, a problem known as thermodynamic inconsistency 
of the approximation. 

Some liquid state theories contain one or more adjustable state-dependent parameters which are tuned 
to achieve thermodynamic consistency between several routes, usually the compressibility and the virial 
ones. Examples include, among other approaches, the modified hypernetted-chain closure [4], the Rogers- 
Young closure [5], the Zerah-Hansen closure [6], the self-consistent Ornstein-Zernike approximation [7], 
the hierarchical reference theory [8], Lee's theory based on the zero-separation theorems [9], the generalized 
mean spherical approximation [10] or the rational- function approximation [11]. 

On the other hand, standard theories do not have fitting parameters and thus they are in general 
thermodynamically inconsistent. An interesting result, however, is that the hypernetted-chain (HNC) 
integral equation provides thermodynamically consistent results through the virial and energy routes, 
regardless of the potential 4>{r) [12]. A similar result has recently been reported [13] in the case of the 
mean spherical approximation (MSA) applied to soft potentials, such as the Gaussian core model 4>(r) = 
eexp[— (r/cr) 2 ]. Therefore, a certain close relationship between the energy and virial routes seems to exist, 
at least for some approximate theories and/or some interaction models. This has been further supported 
by a recent proof on the equivalence between both routes when taking a hard-sphere (HS) limit from the 
square-shoulder (SS) model, regardless of the approximate theory used to describe the structural properties 
of the fluid [14]. 

In a HS liquid, the second term on the right-hand side of Eq. © vanishes, so that the internal energy per 
particle is the same as that of an ideal gas, i.e., it reduces to the kinetic contribution 3/2/3 and is independent 
of density. Moreover, the compressibility factor Z of a HS liquid is independent of temperature. As a 
consequence, the thermodynamic relation © is trivially satisfied as = and it is in principle impossible 
to obtain the EoS of the HS fluid through the energy route. This ill definition of the energy route of a HS 
fluid can be avoided by considering a suitable interaction potential that reduces to the HS one in certain 
limits. The simplest choice for such a potential is perhaps the SS function 



<Pss(r) 



oo, r < a', 



e, a' < r < a, (6) 
0, r > a, 



where e is a positive constant that measures the height of the square shoulder, while the width is given by 
the difference a — a' . The potential © has been studied by several authors in different contexts [15]. For 
the purpose of this paper, it is chosen here because it reduces to the HS interaction potential in several 
special cases. First, in the limit of zero shoulder height, e — > (or, equivalently, in the limit of infinite 
temperature), the potential © becomes that of HS of diameter a'. In the opposite limit of infinite shoulder 
height, e — > oo (or, equivalently, in the limit of zero temperature), one also gets the HS interaction, but this 
time the one corresponding to a diameter a. These two limits are important because one of them is needed 
as a boundary condition when integrating the internal energy over temperature to get Z from Eq. ©. 
Interestingly, the HS potential of diameter a is also recovered in the limit of zero shoulder width, a' — > er, 
at finite e (finite temperature). In the opposite limit a' — > 0, the hard-core part of the interaction has a 
vanishing influence and the SS potential becomes that of so-called penetrable spheres (PS). The PS fluid 
has been extensively studied [16] as an example of bounded potentials describing the effective two-body 
interaction in some colloidal systems [17]. These special limits of the SS potential are sketched in Fig. ^ 
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Apart from that, it is worth recalling that the square-well (SW) potential is described by Eq. (jSJ), except 
that e is negative in that case. From the SW potential one can still get the HS potential of diameter a' 
in the limit |e| — > (or, equivalently, in the limit T — > oo, what then provides the boundary condition for 
the energy route to the EoS) and the HS potential of diameter a in the limit a' — > a, but not either of the 
other two cases {a' — > and |e| — > oo) depicted in Fig. ^ 

As said before, starting from the EoS for SS particles obtained via the energy route, it can be proven [14] 
that when one first takes the limit e — > oo and then the limit a' — > o~, the resulting EoS for HS always 
coincides with the one obtained via the virial route. However, there is no a priori reason to expect that 
this equivalence still holds when the HS limit is taken following a different path. The aim of this paper is 
to clarify this issue and show, by taking a couple of alternative paths, that the energy-route EoS depends 
indeed on the path followed from SS to HS. 

This paper is organized as follows. The virial and energy routes to the thermodynamic properties of 
SS fluids are presented in section [21 The HS limits along different paths are worked out in section The 
results are summarized and discussed in section^ Finally, the paper is closed with the conclusion section. 



2 Virial and energy routes to the equation of state of square-shoulder fluids 

Let us consider a fluid of particles interacting via the pairwise SS potential given by Eq. ©. Therefore, 



e -P<hs(r) 



0, r < a', 

e-P £ , a' < r < a, (7) 

1, r > a. 



We take the "outer" diameter a as fixed and define the relative value A = a' /a of the "inner" diameter 
a'. This quantity < A < 1 parameterizes a family of independent SS potentials. Given a value of A, the 
thermodynamic state is determined by the density and the temperature. In dimensionless units, we can 
measure the density by the parameter r] = (ir/6)pa 3 [18] and the temperature by (5* = f3e = e/ksT = l/T* . 
Moreover, without loss of generality, henceforth the distances will be understood to be measured in units 
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of a. For the SS interaction ([Hjl. the virial and energy equations, Eqs. (Q) and ©, become 



Z&fo, /T; A) = 1 + 47/ A 3 e-^yss(r = A|r?, /?*; A) + (1 - e-"*)i/ss(r = % F\ A) 



(8) 



_3 

2^ 



- + 12r ? e" /3 * J 1 drr 2 yss{r\ V ,/3*;X) 



(9) 



respectively. The notation in Eqs. © and @ makes explicit the dependence of the SS quantities on r/, (3* 
and A. Besides, the superscript v in Eq. (JHJ) is introduced to emphasize that it corresponds to the virial 
route. Making use of Eq. (JSJ), the energy equation © yields the following expression for the compressibility 
factor, 



(r\rj,l3*;X), 



(10) 



where use has already been made of the mapping SS^HS in the infinite-temperature limit (3* — > 0. The 
superscript e affects to Zgs as well as to Zhs since both sides of Eq. (jlOj) must agree in the limit 0* — > 0. 

The virial series expansions for the cavity function and the compressibility factor of the SS fluid are 
defined as 



yssH?7,/3*;A) = l+£>f (r|/?;A)7f, 



:n) 



n=l 



Z S s(r7,/r;A) = l + X>f + i(/?;A>f 



(12) 



n=l 



Inserting the expansion (|12[) into Eqs. (jHJ) and (|10|) . one gets the following expressions for the virial-route 
and energy-route virial coefficients in the SS model: 



bf A) = 4 rA 3 e-^yf_ 2 (A|/3*; A) + (1 - e^>f_ 2 (l|/T ; A)' 



n > 2, 



(13) 



&f ' e (/T; A) = ^X^-D + i2(n - 1) f dft* e~« dr r 2 ^^*. A); n > 2 . 

J0 



(14) 



3 Hard-sphere limits 

Given an (approximate) cavity function yss( r \v> Pi > A), the compressibility factor obtained from Eq. (fTT 
differs in general from the one given by Eq. Q. The question is, does that difference persist in the HS 
limit? To address this point, we need to specify the path followed to get the HS limit. 

Figure El represents the relevant three-dimensional parameter space for SS fluids. The plane a 1 /a = A = 
represents the PS two-dimensional parameter space. The planes fc^T/e = 1/(3* = and a 1 /a = A = 1 
correspond to HS of diameter <r, whose dimensionless properties depend only on the reduced density 
1} = (ir/6)pcr 3 and should be independent of A (on the plane ksT/e = 0) and of temperature (on the plane 
a' I a = 1). The point P represents a SS liquid with a given value of A and at a given thermodynamic state 
(77, /?* ) . From it, one can follow several paths by changing A and/or (3* to reach a HS liquid at the same 
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Figure 2. (Color online) Parameter space for SS fluids. The plane a' /a = represents the PS fluids, while the planes kgT/e = and 
a' /a = 1 correspond to HS fluids. Starting from a given SS fluid (represented by the point P), it is possible to go to the HS fluid at the 
same density by following different paths. In particular, the paths A, A+A', B+B' and C are considered in the paper. 

density 77. Here we will be concerned with paths A (i.e., (5* — > 00), A+A' (/?* — > 00 followed by A — > 1), 
B+B' (A -» followed by (3* -» 00) and C (A -> 1). 



3.1 Path A+A' 

Let us take the zero temperature limit, (3* — > 00 in both sides of Eq. (jlUj) . The result is 



a poo i-l 

Zfain) - ^ s (r?A 3 ) = 12,7—77 J d(3*e~^ J drr 2 y^ 



r\ V ,(3*;X). 



(15) 



This is actually the result obtained through path A, which will be analyzed later on. Before proceeding 
with the subsequent step A', let us rewrite Eq. (|15l) in an equivalent form. First, note that Zhs(??A 3 ) can 
be expanded around 77: 



Z HS (r ] X 3 ) = Z HS (v) + J2 



(-r/) n (l - A 3 ) ra d n 



n=l 



111 



dr\ n 



Zus(r])- 



(16) 



As a consequence, Eq. (|15[) yields 



z mv) + 1 ( "t+i)! 3r i^ z " s(??) = i+ r^r w e-r f dr r2yss(r|?? ' ^ a) ' (i7) 



n=l 



where use has been made of the consistency condition Zhs(0) = 1- Taking now the limit A — > 1 (path A' 
in both sides of Eq. (|17|). and using the physical condition 



lim yss{r\v,P*; A) = yns(r\rf), 

A — > 1 



(18) 



we finally get 



Zbs(v) = l + 47/yHs(%). 



(19) 



This is not but the virial EoS, Eq. (jHJ), particularized to HS. This proves that, no matter which approx- 
imation is used to get 2/ss( r l r 7; P* ', A), the energy EoS coincides with the virial one when the HS limit is 
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reached from the SS fluid following the double path A+A'. The proof presented in Ref. [14] is slightly 
more general since it applies to mixtures and to any dimensionality. 

The equivalence between the energy an virial routes when the path A+A' is followed can also be proven 
at the level of the virial coefficients. We first take the limit j3* — > oo on both sides of Eq. (|14|) (path A) 
with the result 

b n S ' e = J o °° W e ' r [ dr r*yf- 2 {r\P ; A), n > 2. (20) 

Next, the limit A — > 1 (path A') yields 

5 HS,e = b US,v = 4y HS 2(1)> „ > % (21) 



3.2 Path B+B' 

Let us consider now the path B+B' in Fig.[2J which is very different from the path A+A' considered above. 
We will restrict ourselves to the fourth virial coefficient since this is enough to check that the virial and 
energy routes do not coincide now in the HS limit. We start taking the PS limit A —* (path B) on both 
sides of Eqs. (JT3j) and ifrijl. 

^ S ' w C9*)=4(l-e-' J -)^ 2 (l|^), n>2, (22) 



b™> e (P*) = 12(n - 1) / dPl e^* f 1 drr 2 y^_ 2 (r\Pl), n > 2. 
Jo Jo 



(23) 



The exact function U2 S ( r \P*) for the PS model, as well as the corresponding expressions in the Percus- 
Yevick (PY) and HNC approximations, have recently been obtained [19]. In particular, the PY result 
is 



yP- pY (r\n 



;i-e-n : 

35r 



(r - l) 4 (r 3 + 4r 2 - 53r - 162) (^4e _/3 * - l) 6(1 - r) 



+2(r - 2) 2 (r 5 + 4r 4 - 51r 3 - 10r 2 + 479r - 81) (l - e _/3 * J 9(r - 2) 
-(r - 3) 4 (r 3 + 12r 2 + 27r - 6)6(3 - r)] , 

where @(x) is the Heaviside step function. Inserting Eq. (|24l) into Eqs. (|22() and (J23|). one gets 



(24) 



e PY "(r) = |(l-e^ 



35 - 171e" /3 ) , 



(25) 



6 PS-PY, e(/r) = -L (l - e-' 3 *) 4 (907 - 6347e"^) , (26) 

respectively. As expected, the virial and energy routes to the fourth virial coefficient of the PS fluid differ 
in the PY approximation. This difference persists in the HS zero-temperature limit (3* — > oo (path B'), 
namely, 



(27) 
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Figure 3. (Color online) Temperature dependence of the fourth virial coefficient for the PS fluid. The curves correspond to the exact 

result ( — ), the virial and energy routes in the HNC approximation ( ), the compressibility route in the HNC approximation ( ), 

the virial route in the PY approximation (—••—), the compressibility route in the PY approximation (—• — •) and the energy route in 

the PY approximation (■■■). 



tj ,HNC-PY,?j 

However, o 4 



(/?*) = 6 4 " (/?*) [19], in agreement with a general property of the HNC approx- 
imation [12]. Figure El shows the temperature dependence of the exact fourth virial coefficient of the PS 
fluid, as well as the results obtained from the PY and HNC approximations via the virial, compressibility 
and energy routes [19]. 



3.3 Path A 

Equations (|19|) and (|21|) show that the virial and energy routes are always equivalent when the HS limit is 
taken through the double path A+A'. On the other hand, Eq. l|27j) shows that this equivalence is generally 
broken when the chosen path is B+B'. While in the path A+A' one ends with A = 1, in the path B+B' 
the first step is A — > 0. Therefore, it might be reasonably expected that if one directly goes from SS to 
HS through the path A (see Fig. EJ) the resulting energy-route EoS for HS artificially depends on A. To 
illustrate this, let us consider the following toy approximation for the function y^{r\f3*; A): 



SS-toy 
2/2 



(r|/T;A) 



y P 2 S ~ PY (r\n 



+ y r- PY (r/X)X e 



v-2 



PS - pY (r|/3*)y 2 HS - pY (r/A)A 6 



y HS-PY (r) 



(28) 



where y PS_PY (r|/3*) is given by Eq. (glj) and yf Y (r) = lim y g»_» Q0 y PS (r|/3*) 
(|28jl reduces to the PY results in the four limits indicated in Fig. namely 



The toy approximation 



^lim % SS - toy (r|/3*;A) 



, 2 HS " PY (r/A)A 6 , 



(29) 



lim y 2 SS - toy (r|/T; A) = km yf- toy (r\P* ■ A) = y 2 HS " PY (r), 

A — ^1 



(30) 



(31) 
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Figure 4. (Color online) Plot of the HS fourth virial coefficient as a function of A = o~'/o~. This coefficient is obtained by starting from 
the energy-route coefficient for the SS fluid in the approximation 1281 and then taking the zero-temperature limit (path A). 

Prom that point of view, y2 S_toy (r|/?*; A) can be seen as a simplified version of the true yf s " PY (r|/?*; A) 
provided by the PY approximation. In any case, the point here is not how accurate or how close to the 
PY function the toy approximation (|28[) is, but to illustrate the sensitivity of the energy route to the fixed 
value of A. 

According to Eq. (|2()j) . the energy-route fourth virial coefficient when the path A is followed becomes 



,HS-toy,e _ 
°4 — 



36 

r^A9 



dr r z 



d/3* e-P*yf- toy (r\p*;\). 



(32) 



The integration over (3* is straightforward, while in the integration over r one needs to distinguish the 
cases 0<A<|,^<A<^ and 1 < A < 1. The artificial dependence of ^ s_toy,e on A is shown in Fig. 0J 
The extreme points of the curve agree with the results previously obtained, i.e., 6^ S_PY ' e = b^ S ~ PY,v = 16 
along the path A+A' but 6^ S_PY ' 6 = 1814/175 ^ b^ S ~ PY ' v along the path B+B'. It cannot be ascertained 
at this point whether the non-monotonic behavior in Fig. |1] is an artifact of the toy approximation (|28|) 
or is a feature also shared by the PY approximation. In any case, the important issue here is that the 
equivalence between the energy and virial routes is only reached, in general, if A — > 1. 



3.4 Path C 

As Fig. [21 illustrates, the HS fluid can also be reached from the PS fluid by keeping the temperature 
constant but shrinking the shoulder (path C). However, this does not provide any information about the 
energy-route EoS of the HS system. Taking the limit A — > 1 on both sides of Eq. (|ll)jl. one simply gets 
Z|j S (r/) = Zy S (r/). In any case, an interesting equation is obtained by first differentiating with respect to 
A and then taking the limit A — > 1. The result is 



dZl s ( V ,P*;X) 



OX 



3?? 



A=l 



d_ 

dr] 



-13* 



Z v us(v) 



(33) 



4 Summary and discussion 



The aim of this paper has been to investigate the possibility of circumventing the ill definition of the 
energy route to thermodynamics for HS fluids by first considering SS fluids and then taking the adequate 
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limits. The SS interaction potential is particularly appropriate because it is simple and yet reduces to the 
HS model in several independent limits (see Fig. 1). The high-temperature limit (T* — > oo or, equivalently, 
e — > 0) is important since it provides the necessary boundary condition to get the compressibility factor 
by integrating the internal energy over temperature [see Eqs. (jSJ, l)lflj> an d p4|>]. 

Imagine that an approximate cavity function j/ss( r I 7 ?; P* ', A) for the SS liquid is known (either analyti- 
cally or numerically). Then, Eqs. Ijlfljl and/or (j!4j) can be used to assign a meaning to the energy-route 
compressibility factor for HS, Z^ s (r/), or to the associated virial coefficients, b n ' e . In order to do so, the 
limit of vanishing shoulder width at finite temperature (path C in Fig. EJ) is useless. On the other hand, 
well defined results are obtained by taking the zero-temperature limit (3* — > oo (path A), as shown by 
Eqs. (|15|) and ()2()ll . Although the above limit is taken at fixed width A, the final result should, on physical 
grounds, be independent of A. However, the approximate nature of yss(f\Vi A) gives rise, in general, to 
an inconsistent dependence of Z^ s (rj) on A, as illustrated by Fig. 0]in the case of the toy approximation 

J2EJ). 

The artificial dependence of Z^ s (rj) on A suggests to further take the limit A — > 1 as the most "sensible" 
path to go from SS to HS (path A+A'). In this way, the HS fluid is reached twice, first by decreasing the 
temperature (or, equivalently, increasing the shoulder height, e — * oo) and then by shrinking the shoulder 
width. As proven here and in Ref. [14] , the resulting EoS coincides exactly with the one obtained directly 
from the virial route, i.e. Z^ s (rj) = Z^ s (rj), no matter which approximate theory is used. 

As said above, the anomalous A-dependence of Z^ s (rj) when the path A is followed has been illustrated by 
considering a simple toy approximation, Eq. (j2SJ), and thus it might be conjectured that such a dependence 
disappears when a most "respectful" theory is taken into account. However, this is not the case, at least 
for the PY theory. When the zero-temperature limit (3* — > oo is taken at A = (i.e., from the PS model, 
path B+B'), the resulting value of the fourth virial coefficient (b^ " PY,e ~ 10.37) strongly differs from the 
value obtained in the opposite limit A — > 1 (6^ S_PY ' e = b^ S ~ PY ' v = 16). 

Of course, all the routes to thermodynamics self-consistently agree if the exact function yss( r l r ? ! P*\ A) 
is considered. Moreover, the HNC theory is known to yield consistent thermodynamic properties via the 
energy and virial routes, for any interaction potential [12] . This is explicitly verified in the case of the fourth 
virial coefficient for PS fluids [19], as shown in Fig. |31 Therefore, a plot (similar to that of Fig. 0} of the 
fourth virial coefficient obtained from path A would show the constant value 6 4 ,e = 6 4 ' = 28.5. 
The interesting question is, which is the necessary and sufficient condition to get thermodynamic properties 
independent of A when going from SS to HS through path A? To address this question, define the quantity 



poo rl 

*asfa;A) = 12»7 / d/Te"^ / dr 
Jo Jx 



r 2 y ss (r\ V ,t3*;X). (34) 



Then, either from Eq. (|15|) or from Eq. Q using the thermodynamic relation u = d((3f)/d(3, where / is 
the free energy per particle, one gets 

<f> ss ( V ;\) = <p e HS (r l )-<p e m (r l \ s ), (35) 

where (p = f3f ex , f ex being the excess free energy per particle. Therefore, the sought necessary and sufficient 
condition is that the quantity &ss(v'i A) must be equal to a function that only depends on rj plus a function 
that only depends on ry' = ?7A 3 . The fact that both functions are actually the same, except for a sign, is a 
consequence of the trivial property $ss( r ?i 1) = 0. In differential form, the condition (|35|) becomes 

|. A |. $SS( „ ;A) = 3^|-* SS (,;A). (36) 

If this condition is fulfilled, then Eq. (|33|) implies that ^fjg(^) = < I ) ss( 7 ?; 0). The independence of yfisl 7 ?) 011 
A implies that the energy and virial routes become equivalent, as discussed in the text. This can be easily 
checked from Eq. Q35JI by differentiating both sides with respect to A and then setting A = 1. Taking into 
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account the thermodynamic relation Z = pd((3f)/dp, one then gets Z^ s (r]) — 1 = Z^ s (r]) — 1 = ^rjyi{s(l\r]). 
5 Conclusion 

The question posed in the title of the paper is only meaningful if the energy route is understood by starting 
from an interaction potential (p(r) which encompasses the HS model in certain limits. Since, according to 
Eq. ©, the obtention of the compressibility factor from the internal energy requires an integration over 
temperature, it is necessary that the potential (fi(r) becomes equivalent to that of HS (or negligible) in the 
limit T — > oo. Next, in order to get a non trivial result, 4>(r) must also become indistinguishable from HS 
in an independent limit, for instance T — > 0. The results presented in this paper by taking cp(r) = (f>ss( r ) 
suggest that, in general, either the energy-route EoS for HS fluids is identical to the virial-route EoS or 
the former is not unique but depends on the path followed to reach the HS fluid. 

A final comment is in order. When considering the energy route in the case of SW fluids, it is usual 
to fix the boundary condition at T — > oo by freely choosing a convenient form for Zhs(ji), such as the 
Carnahan-Starling EoS [3]. However, this must be done with caution since the resulting EoS would become 
inconsistent when making the change e <-> — e, in which case </>sw( r ) <^ss( r )- From that point of view, it 
would be more consistent to take the virial form Z^ s (j]) corresponding to the approximation (PY, HNC, 
MSA, . . . ) being used to describe the SW fluid. 

This work has been supported by the Ministerio de Education y Ciencia (Spain) through Grant No. 
FIS2004-01399 (partially financed by FEDER funds) and by the European Community's Human Potential 
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